Measuring the kernel of time-dependent density functional theory with X-ray 
absorption spectroscopy of 3d transition metals 
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The 2p — 3d core hole interaction in the Z/2,3 absorption spectra of the 3d transition metals 
is treated within time-dependent density functional theory (TDDFT). A simple three- level model 
explains the origin of the strong deviations from the one-particle branching ratio and yields matrix 
elements of the unknown exchange-correlation (XC) kernel directly from experiment. 
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Ground-state density functional theory (DFT) is well- 
established for atoms, molecules, and solids l\. But 
ground-state DFT produces only a one-particle picture 
of the electronic transitions in matter, neglecting inter- 
actions between excitations. In the optical regime, time- 
dependent DFT (TDDFT) has enjoyed recent success 
M 0] m describing such dynamic exchange-correlation 
(XC) effects. The spectroscopic properties of matter in 
the X-ray regime are substantially governed by dynami- 
cal many-body effects involving the creation of a localized 
core hole 0,0,0,110. While GW calculations and the 
Bethe-Salpeter equation can be used 0], these are com- 
putationally demanding. The less expensive methodol- 
ogy of TDDFT is now being developed for these effects 



We analyze this approach to the X-ray absorption of 
itinerant systems like the £2,3 absorption of 3d transition 
metals (TMs), i.e., exciting a photoelectron from the lo- 
calized 2p core states into the 3d band. L2.3 X-ray ab- 
sorption spectra (XAS), especially of early 3d TMs, suf- 
fer from core-hole correlation effects 0. Schwitalla and 
Ebert Jjj] applied TDDFT linear response theory to calcu- 
late the XAS of the 3d TMs. Using a local approximation 
to the frequency-dependent XC kernel, as proposed by 
Gross and Kohn [ljj, they qualitatively reproduced the 
trend of the branching ratios. However, whenever DFT 
is applied in a new regime, a difficult question arises: Are 
the existing functional approximations sufficiently accu- 
rate in this new regime? And how does one separate XC 
errors from those due to the practical approximations 
needed for realistic calculations? In Ref. |8j, the limita- 
tions of the Gross-Kohn approximation for this problem 
are shown, and a new approximation suggested. But the 
true value of DFT is in constructing one XC approxi- 
mation that covers many situations, in order to build-in 
knowledge of the underlying physics. Other TDDFT ap- 
proximations, such as Vignale-Kohn would need to 
be inserted into their codes in order to be tested. 

Our approach here is different, and is based on the phi- 



losophy of Ref. 12] . That work examined the TDDFT re- 
sponse when excitations are not strongly coupled to each 
other. A useful series was developed in the strength of the 
off-diagonal matrix elements, relative to the frequency 
shifts induced by diagonal terms. The leading term yields 
the single-pole approximation (SPA) 2] , which has proven 
very useful in understanding TDDFT corrections to the 
one-particle picture. It even yields an immediate esti- 
mate of the XC kernel, but only if excitations are well- 
separated, a criterion rarely realized in practice fl3j | . 

However, the same philosophy applies to cases of two 
levels strongly coupled to one another, but weakly cou- 
pled to the rest of the spectrum. We call this the three- 
level or double-pole approximation (DPA), cf. Fig ^ 
Moreover, the £2,3 absorption of 3d TMs provides an 
ideal example of two transitions much closer to each other 
than the rest of the spectrum. With this in mind, we 
experimentally measured the branching ratios and level 
splittings of the 2p 3 / 2 (^3) and 2^1/2(^2) core states, and 
now deduce off-diagonal matrix elements of the unknown 
XC kernel. Because we can also compare with the one- 
particle Kohn-Sham (KS) spectrum, we can also deduce 
the diagonal matrix elements. The large deviation of 
branching ratios from their single-particle values is due 
entirely to the effect of core-hole interaction on spin-orbit 
coupling. Thus the DPA to TDDFT explains the ob- 
served shifts and oscillator strengths, and also provide 
benchmarks for future XC kernel approximations. We 
believe this is the first experimental measurement of a 
matrix element of the XC kernel of TDDFT. 

Consider a system of electrons subject to a small 
frequency-dependent perturbation. By virtue of TDDFT 
the corresponding linear density-density response func- 
tion x is related to the response function Xs of non- 
interacting particles via the Dyson-type equation |lfll | 

X(r,r',w) = Xs(r,r',uj) 

+ J d 3 xj dV Xs (r,x,w)X(x,x>) X (x>>).(l) 



2 



The kernel K(r,r',ui) consists of the bare Coulomb 
interaction and the frequency-dependent XC kernel 
/ xc (r,r',w): 



E 3 (3d) 



tf(r,r',u,) = — 6 +f xc (r,r',u). 
r — r 1 



(2) 



The exact XC kernel describes, among other many-body 
effects, the core-hole interaction with the photoelectron. 
Neglecting K, the spectrum would reduce to the bare KS 
single particle spectrum represented by \s- In XAS, the 
deviations produced by K are called core-hole correlation 
effects. The response function Xs is given in terms of the 
ground-state KS orbitals tpj (spin-saturated) by 

X s {r,r',Lj) = 2> (rij-n k ) J - — (3) 

v - ujj k + it] 

where rij,nk denote the Fermi-occupation factors and 
u)jk are the KS orbital-energy differences. For a single- 
particle transition q (q = k — > j) define $ 9 (r) := 
(pk(r)(p1j(r). The exact density-response function \ nas 
poles at the true, correlated, excitation energies Qj, 
which can be found by solving [T3 |: 



9' 



1,3 1 



(4) 



where the matrix is 

.2 



wv(n) = uj; s qq , +A^cj q Lj' q K qqf (ii), 

K qq ,{n) = Jd 3 r f d 3 r'$ q (r)K(rr'Q,)$ q ,(r'). (5) 



The eigenvectors yield the oscillator strengths [lj] via 



\x T S-^ 



(6) 



where S qq ? — 5 qq i ^ujq~ and Xq IS 3j column of the KS 
dipole matrix elements. This eigenvalue problem rigor- 
ously determines the excitation spectrum of the interact- 
ing system, but the quality of the results in any practical 
calculation depends crucially on the approximation em- 
ployed for the XC kernel. 

In the £2,3 XAS of the 3d TMs, the description of 
the electron core-hole interaction may be simplified by 
the assumption that the relativistic spin-orbit coupling 
(SOC) in the 3d band states (~ 0.05 eV) is small com- 
pared to that of the core states (several eV) and can be 
neglected. This means that the oscillator strengths fj of 
these levels are all about equal, as their KS orbitals are 
essentially identical. Since, in this limit, the absorption 
area is proportional to the oscillator strength, weighted 
statistically according to the manifold of the j = 3/2 
and j = 1/2 subshells, the branching ratio of the KS sys- 
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tern is B KS = A 3 / 2 /(A 3 / 2 + A 



h/2) 



2/3, where Aj is 



FIG. 1: Schematic illustration of the DPA model. The model 
describes the shifts of the excitation energies (uncorrelated u>i 
and correlated fl j) and the changes in corresponding oscillator 
strengths /j in the presence of an excited core hole. 



the area under the peak of the j-th subshell. (^46 initio 
calculations of the L 2: 3 XAS without core-hole correla- 
tion effects based on the fully-relativistic spin-polarized 
Korringa-Kohn-Rostoker band structure formalism yield 
branching ratios very close to £?ks [H[3.) Here we 
replace all dipole-allowed transitions cjjk from a particu- 
lar absorption edge into the 3c? band by a single particle 
transition, as illustrated in Fig. 2] 

In Fig. |2 we show our experimental isotropic XAS 
for the 3d TM with almost empty 3c? bands taken from 
Fe/TM/Fe sandwiches with TM = Ti,V, Cr and bulk- 
like Fe. The data were recorded at the UE56-1/PGM 
beamline at BESSY (for details, see Ref. 13). The edge 
jumps are normalized to unity. From these spectra and 
their absolute energy dependence, the excitation energies 
Cl q -i at the L3 edge and fl q = 2 and the L 2 edge are de- 
termined. For the quantitative analysis of the branching 
ratio B, we very carefully determined the L 2 3 absorption 
areas Aj . This determination has the advantage that B 
becomes independent of the different L3 and L 2 lifetime 
broadening and experimental resolution. Note, that the 
proper experimental intensity is given by the area and 
not by the height of the resonance. To determine the 
correct area of the L3 and L 2 resonances the continuum 
contribution is removed (e.g. gray line for Fe in Fig. [21 
H3)- Since the 2p SOC decreases towards lower atomic 
numbers the deconvolution is more complicated for the 
early 3d TMs Ti, V, and Cr because of the strong L 2j 3 
overlap. The areas have been fitted using the Fe absorp- 
tion spectrum as a background simulation underneath 
the L 2 edge. This fit appears justified, since the L 2 on- 
sets for the early 3c? TMs follow systematically the energy 
dependence of the L3 edge of Fe. The experimental re- 
sults are set out in Table HJ The energy resolution of the 
experimental spectra (Fig. |2J) is approximately 0.5 eV. 
Consequently, the widths of a few eV in the L3 and L 2 
resonances are mostly due to lifetime effects. 

In the case of Fe the L 2 absorption is approximately 
half of the L3 peak, in agreement with the KS prediction. 
However, the branching ratios for the other 3c? elements 
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FIG. 2: The experimental isotropic absorption spectra (solid 
line) at the Z/2,3 edges are shown for the early 3d TMs Ti, V, 
and Cr versus Fe. The edge jumps are normalized to unity 
for direct comparison. The continuum in the experimental 
spectrum is simulated by a two-step function as shown for 
Fe (dashed-dotted line). The treatment of the core hole red- 
shifts the independent particle spectrum (dotted line) and 
changes the statistical branching ratio in the correlated spec- 
trum (dashed line) as revealed by the DPA model. 



differ significantly from this. In particular, Ti has an L 2 
peak that is even larger than its L 3 absorption. Thus, the 
experimental branching ratios cannot be interpreted in 
terms of KS orbitals, suggesting strong electron core-hole 
interactions. In the language of TDDFT, there must be 
significant off-diagonal matrix elements in Casida's equa- 
tions, describing the influence of the electron core-hole 
interaction on the L2.3 XAS. (If only diagonal elements 
are considered in Eq. (0J, the eigenvalues are shifted 
but the eigenvectors are not rotated, and the oscillator 
strengths retain their KS values |l2j.) 

However, a fully numerical solution of the equations 
is not needed, as we know there are only two dominant 
transitions, so the electron core-hole interaction can be 
analyzed within the DPA model. For only two transi- 
tions, we solve Eq. (0J exactly: 



^ sin(c?/2)) 2 
IZi sin(c?/2)) 2 



(7) 



fi = (V7mCos(0/2) 
H = {yfhfl cos(#/2) 
for the interacting oscillator strengths fi, where 

taae = 2Wis/{Wm-W n ) (8) 
is a mixing angle that represents the strength of the cou- 



pling between the two transitions. This yields: 

„ 4- v // s , 1 / s . 2 qiq 2 + (/s.2 - /s,i)(<?2 - ?i) /n s 
cos 9 = — — — — — - — - — r (9) 

(/s,i + fs,a)[qi + Q2) 

where / Sii denotes a Kohn-Sham oscillator strength and 
qi = Bp2, Q2 = (1 — B)pi, where pi the multiplicity 
of the initial state of transition i. In our case, setting 
/s,i = fs,2 = fo, Eq. simplifies to 



h = f (l -sin9), h = h{l 



sin ( 



(10) 



and since the absorption depends only on the oscillator 
strengths, weighted statistically, Eq. © yields 



sin 9= (2 -SB)/ (2- B). 



(11) 



The branching ratio directly determines the mixing an- 
gle! We can even extract directly the off-diagonal matrix 
element of the Hartree-XC kernel itself. Inserting the 
matrix elements into Eq. (JSJ, and neglecting all small 
differences (e.g., between KS and exact transition fre- 
quencies) which are a few eV in several hundred, 



K 12 = sin 9 Aft/4 . 



(12) 



Thus knowledge of both the branching ratio and the level 
splitting AQ = O2 — S^i yields an experimental determi- 
nation of the off-diagonal matrix element of the XC ker- 
nel. Lastly, given the ground-state KS energy levels, we 
can also recover the diagonal elements: 

K j:i = (ft 1 +ft 2 )/4+(-l) J 'Aft cos(9 /4-w,-/2 j = 1,2. 

(13) 

These results are listed in Tabled the corresponding the- 
oretical DPA spectra are presented in Fig. [21 

This analysis yields a simple interpretation of the ob- 
served spectra. First, imagine there was no SOC. Then 
there would be a single p-level, and SPA applies. The 
diagonal matrix element of / HX c is the well-known core- 
hole interaction that shifts the transition frequency from 
its KS value. In the presence of spin-orbit splitting, both 
levels are shifted by similar amounts, about 5-7 eV. 

Much more importantly, a new effect occurs, which is 
that the core-hole interaction couples the two KS tran- 
sitions together, altering their branching ratio. The 
much smaller off-diagonal core-hole interaction (about 
1/2 eV) produces the large deviation from the single- 
particle branching ratio. Although the matrix elements 
are about 5 times smaller than their diagonal counter- 
parts, in Ti they reverse the relative sizes of the peaks. 
This effect can be thought of as simple level (or in this 
case, transition) repulsion, as the two transitions near 
one another. Eq. (|12f> shows that the true measure of 
coupling is 4|ifi2|/Af2, which is growing from Fe to Ti 
only because the 2p SOC Aft is shrinking, not because of 
increased interaction. We also note that SPA |2| (which 
can be recovered in all results by setting 9 = 0) yields 
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TABLE I: Excitation energies in eV obtained from KS 
calculations (wf" S ) and from experiment (fij), experimental 
branching ratio B and matrix elements Kij. The experimen- 
tal error of fi; is below 10 -3 , the one of B in the order of 
1 %. 



3d TM 


, , KS 


, , KS 


Qi 


n 2 


B 


K u 


K 22 


K 12 


22 Ti 


460.8 


467.5 


455.4 


461.0 


0.47 


-2.57 


-3.34 


0.54 


23 V 


519.1 


527.7 


513.6 


520.4 


0.51 


-2.65 


-3.73 


0.54 


24 Cr 


580.3 


590.3 


575.1 


583.6 


0.56 


-2.55 


-3.40 


0.47 


26 Fe 


711.3 


724.6 


706.7 


719.5 


0.70 


-2.29 


-2.55 


-0.25 



2Kjj — Qj — ujj and is highly accurate for the diagonal 
elements. Thus the shifts are simply interpreted as di- 
agonals of K, while the branching ratios are a sensitive 
determinant of off-diagonal elements. Lastly, for very 
small splitting (AS1 <§; 4.K12), s'md — > 1 and, from Eq. 
O, all weight goes into the L2 peak. From Eas. (|12ll3fl : 

AO = 2yJ (Auj/2 + K 22 - K n ) 2 + 4Kf 2 . (14) 

Thus 4| 12 1 is the minimum level splitting, and occurs 
with the L2 peak much larger than L3. 

The success of DPA shows that very little effort beyond 
a ground-state DFT calculation is needed to compute 
these spectra in TDDFT. One only needs to integrate a 
given approximation to the XC kernel for the two diago- 
nal matrix elements, and one off-diagonal. Furthermore, 
applying our analysis in reverse to the calculated ALDA 
and RPA results of Fig. 1 of Ref. || , we find that the 
success of the suggested approximation (which implies 
using ALDA for diagonal elements, and neglecting XC 
for the off-diagonal elements, i.e., RPA) also implies that 
ALDA is accurate for the peak shifts alone, while RPA is 
accurate for the branching ratio predicted by Eq. (|12[) . 
once the experimental shift is used. 

Why is DPA justified for these systems? There are two 
clear sources of error. On the one hand, there are transi- 
tions to other levels to consider, but these are well-known 
to have small effects For example, the transitions 

to 4s have much smaller oscillator strengths, while other 
transitions are included in the background, which has 
been subtracted. Of greater concern might be the lack of 
resolution of the (slightly) split 3d levels, which yield a 
9x9 Casida matrix problem of allowed transitions. How- 
ever, just as DPA reduces to the single-pole approxima- 
tion of Ref. when levels are too close to be resolved 
[l9) . we expect the full solution of the 9x9 problem to 
collapse to the DPA results when the individual d-levels 
cannot be resolved. This will only be true for the early 
TM's, in which most of the 3d-levels are unoccupied. 

In summary, we have used TDDFT to understand the 
XAS of 3d transition metals by deriving a double-pole 



approximation. The main features observed in the ex- 
periments can easily be explained by assuming that the 
spectrum is dominated by two strongly coupled poles via 
the 2p — 3d core hole interaction. This shows that, for 
the beginning of the 3d series, the reduced 2p-SOC is 
responsible for the strong variation of the branching ra- 
tio, not strong interactions between the transitions. Our 
analysis does not replace a full TDDFT calculation of X- 
ray absorption spectra. Rather, for the very specific case 
of spectral regions dominated by two poles it provides, 
on the one hand, a transparent picture of the changes of 
spectral weights in particular for the early 3d TMs, and 
on the other, a straight-forward route to testing approx- 
imate XC kernels against experimental data. 
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